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In this paper, a modified Eddington-inspired-Born-Infeld (EiBI) theory with a pure trace term 
OiiuR being added to the determinantal action is analysed from a cosmological point of view. It 
corresponds to the most general action constructed from a rank two tensor that contains up to first 
order terms in curvature. This term can equally be seen as a conformal factor multiplying the metric 
This very interesting type of amendment has not been considered within the Palatini formalism 
despite the large amount of works on the Born-Infeld-inspired theory of gravity. This model can 
provide smooth bouncing solutions which were not allowed in the EiBI model for the same EiBI 
coupling. Most interestingly, for a radiation filled universe there are some regions of the parameter 
space that can naturally lead to a de Sitter inflationary stage without the need of any exotic matter 
field. Finally, in this model we discover a new type of cosmic “quasi-sudden” singularity, where the 
cosmic time derivative of the Hubble rate becomes very large but finite at a finite cosmic time. 

PACS numbers: 98.80.-k, 98.80.Jk, 04.50.Kd, 04.20.Dw 

Keywords: modified theories of gravity, early time universe, Eddington-inspired-Born-Infeld theories, cosmic 
singularities 


I. INTRODUCTION 

Undeniably, Einstein’s theory of general relativity 
(GR) has been an extremely successful theory for around 
a century [l|. However, the theory is expected to break 
down at some points at very high energies where quantum 
effects are expected to become crucial, such as in the past 
expansion of the Universe where GR predicts a big bang 
singularity Q. This is one of the motivations for look¬ 
ing for possible modified theories of gravity, which are 
hoped to not only be able to preserve the huge achieve¬ 
ments of GR, but also to shed some light on smoothing 
the singularities predicted in GR. Such theories could be 
seen as effective/phenomenological approaches of a more 
fundamental quantum theory of gravity. 

In 1934 Born and Infeld proposed a non-linear action 
for classical electrodynamics, which is characterised by 
its grand success in solving the divergence of the self¬ 
energy of point-like charges Q. This action for electro¬ 
magnetism has an elegant determinantal structure which 
reads: 



with and being the metric tensor and the electro¬ 
magnetic field strength. Here k is a dimensional constant. 
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Note that the Born-Infeld action (HH) recovers Maxwell 
action for small amplitudes. 

Since the proposal of the Born-Infeld action, modified 
theories of gravity with a Born-Infeld-inspired action, ini¬ 
tiated by the pioneering work of Deser and Gibbons 0 , 
have received much attention; c.f. Refs. [HI 3 . These 
theories of gravity not only maintain the properties of 
GR for small curvatures, but provide various interesting 
deviations from GR at high curvature regimes. Most of 
these works start with a general gravitational action of 
the form 


^det — d \9^,u U (R; Ra/3 ; R/ii/Q/3) | j (^■^) 


where represents a linear Ricci term plus higher 
order curvature terms containing R, R^,j, and R^uap- 
Rather than amending the higher curvature terms in 
we can also modify the action by multiplying the metric 
term in the determinant by a factor 1 -|- aiR, where ai is 
a constant with length square dimensions. This approach 
does not lead to a violation of the requirements in the low 
curvature limit, i.e. we can still recover GR at low ener¬ 
gies. Actually, a gravitational action with this pure trace 
term has been considered in Ref. Q within a pure metric 
formalism. The theory thus inevitably suffers from the 
presence of troublesome fourth order derivatives in the 
field equations or ghost instabilities Q. 

In order to keep the theory free from aforementioned 
problems, alternative theories formulated within the 
Palatini formalism and teleparallel representation have 
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been widely studied in Refs. For example, a 

theory constructed upon the Palatini approach, which 
is dubbed Eddington-inspired-Born-Infeld theory (EiBI) 
(see Ref. M), has recently attracted a lot of attention 
and has been studied from both astrophysical and cosmo¬ 
logical points of view [ 1 M 3 . The EiBI theory is shown 
to be able to cure the big bang singularity for a radia¬ 
tion dominated universe through a loitering effec10 or a 
bouncc0 in the past, with the coupling constant k being 
positive or negative, respectively. The ability of the the¬ 
ory to smooth other cosmological singularities in a phan- 


tom dominated universe was also studied in Refs. [If 
Interestingly, in Refs, the authors showed that the 

bouncing solutions for negative k are robust against the 
changes of the Lagrangian through an additional f(R) 
term or some functional extensions (see as well Ref. |5l| 
for another generalized grav itational theory related to 
massive gravity and Ref. (52j| for the tensorial perturba¬ 
tions of a further generalized gravitational theory within 
the Palatini formalism.). However, it should be stressed 
that the amendments through the addition of a pure trace 
term to the determinantal action have never been consid¬ 
ered so far. Besides, the EiBI theory with a negative cou¬ 
pling constant k was also shown to suffer from instab ility 
problems due to the imaginary effective sound speed . 

On the other hand, a recently proposed determinan¬ 
tal gravity formulated within the teleparallel represen¬ 
tation was shown to be able to cure the big bang sin¬ 
gularity in the past evolution of the Universe through a 
de Sitter inflationary phase [l^. Considering the widest 
generalization, the author added a pure trace term into 
the Lagrangian of the form of QfiuT, where T is the 
Weitzenbock invariant [s^. In our previous work, we 
exhibited that this theory contains cosmological singu¬ 
larities for some parameters of the model, including the 
emergence of some cosmological singularities from purely 
geometrical effects (without the need of exotic matter) 


As far as we know, the gravitational actions with a 
pure trace term added to the Born-Infeld determinantal 
structure have never been considered within the Palatini 
approach in the literature. Eurthermore, we expect the 
emergence of interesting cosmological solutions with the 
addition of a pure trace term because it is expected in the 
teleparallel version. Based on these motivations, in this 
work we will consider a modified EiBI theory with a pure 
trace term added to the determinantal action, and anal¬ 
yse its cosmological implications. Eor simplicity, in this 
work we will assume a homogeneous and isotropic uni¬ 
verse filled with a perfect fluid with a constant equation 


^ The Universe starts its evolution with a minimum size and at an 
infinite cosmic time in the past, before it enters the standard GR 
expansion [T3| . 

^ A bouncing universe is a universe with a minimum or a maximum 
scale factor such that after a contracting phase, an expanding 
phase happens or the other way around. In this kind of model, 
the big bang is substituted by a bounce. 


of state. Because the field equations are complicated, we 
will follow a method similar to that used in Ref. [^ to 
demonstrate the results graphically. 

In this paper, we will follow Ref. to character¬ 

ize the cosmological singularities by the behavior of the 
Hubble rate and its cosmic time derivative at the singular 
event: 0 

• A big rip singularity takes place at a finite cosmic 
time with an infinite scale factor, where the Hubble 
parameter and its cosmic time derivative diverge 

• A sudden singularity takes place at a finite cosmic 
time with a finite scale factor, where the Hubble 
parameter remains finite but its cosmic time deriva¬ 
tive diverges [H, HI, • 

• A big freeze singularity takes place at a finite cos¬ 
mic time with a finite scale factor, where the Hub¬ 
ble parameter and its cosmic time derivative di¬ 
verge [HI, [fiSl - lTlj l . 

• A type IV singularity takes place at a finite cosmic 

time with a finite scale factor, where the Hubble 
parameter and its cosmic time derivative remain fi¬ 
nite, but higher cosmic time derivatives of the Hub¬ 
ble parameter still diverge [HI, HI, . 

• A little rip event takes place at an infinite cosmic 
time with an infinite scale factor, where the Hubble 
rate and its cosmic time derivative diverge Hlii- 

[73. 


• A little sibling of the big rip takes place at an 
infinite cosmic time with an infinite scale factor, 
where the Hubble rate di verg es, but its cosmic time 
derivative remains finite [80j. 

Our results are clearly shown in Table U where we com¬ 
pare them with the original EiBI model [l3- As can 
be seen the model we are proposing can provide smooth 
bouncing solutions which were not allowed in the EiBI 
model for the same EiBI coupling (k > 0). Most in¬ 
terestingly, for a radiation filled universe there are some 
regions of the parameter space that can naturally lead 
to a de Sitter inflationary stage without the need of any 
exotic matter field. Einally, in this model we discover a 
new type of cosmic “quasi-sudden” singularity, where the 
cosmic time derivative of the Hubble rate becomes very 
large but finite at a finite cosmic time. 

This paper is outlined as follows. In section H, we 
briefly introduce the basis of the modified EiBI theory 
with the addition of a pure trace term, including its ac¬ 
tion, field equations, and the low curvature limits of the 


® For an alternative classification of cosmological singularities see 
Refs. 
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theory. In section III, we assume a homogeneous and 
isotropic universe filled with a perfect fluid with a con¬ 
stant equation of state, then follow a similar approach to 
that used in Ref. to derive a parametric Friedmann 
equation. In section IV, we exhibit the evolution of the 
Universe by graphically showing the Hubble rate as a 
function of the energy density under different assump¬ 
tions of the parameters characterising the theory. To 
analyse the evolution of the Universe at the very early 
time, we then confine ourselves to a radiation dominated 
universe in our analysis of the modified EiBI theory. We 
finally present our conclusions in section V. 


II. PROPOSED MODEL: ACTION AND FIELD 
EQUATIONS 


As mentioned in the introduction, in this paper we will 
add a pure trace term, which takes the form of gp,uR, 
to the EiBI determinantal Lagrangian. This term can 
equally be seen as a conformal factor multiplying the 
metric Therefore, the action of this generalized EiBI 
theory is 



K 




\g^Lv + 




( 2 . 1 ) 


where 


= aRf^^{T) + Pg ( 2 . 2 ) 


The theory is formulated within the Palatini formalism, 
in which the metric and the connection T are treated 
as independent variables. In addition, Rfj.u{T) is chosen 
to be the symmetric part of the Ricci tensor and the con¬ 
nection is also assumed to be torsionless. Note that g is 
the determinant of the metric and stands for the mat¬ 
ter Lagrangian, where matter is assumed to be coupled 
covariantly to the metric g only, a, /3 and A are dimen¬ 
sionless constants. The parameter k is a constant with 
inverse dimensions to that of a cosmological constant (in 
this paper, we will work with Planck units SttG = 1 and 
set the speed of light to c = 1). 

In the low energy limit (k —>■ 0), the gravitational ac¬ 
tion inD becomes 




- / ct^XyZ—g {a + ‘if3)R — 2A 


-f (^Q!^ -I- 0/3 -|- ‘ip’^^K.R? 


(2.3) 


where the effective cosmological constant is defined by 
A = (A — 1 )/k. Therefore, the dimensionless constants a 
and /3 should satisfy: 


EiBI theory is regained when 0=1 and /3 = 0. On the 
other hand, this theory becomes an R? theory with its 
gravitational action part being 

Sg\^=o = IJ d^xy/^(^R+^R^ - 2A), (2.5) 

when 0 = 0 and /3 = 1/4 [8l|. Note that this constitutes 
the sole of /(i?)-like action that can be derived from the 
determinantal structure, which is also valid for /(T)-like 
action in the teleparallel representation [T^. Therefore, 
the dimensionless constants o and j3 in this theory can be 
used to quantify the extent of the interpolation between 
the Palatini RZ theory and the EiBI theory. 

Actually, one can also add the so called zeroth or¬ 
der curvature term; i.e., 'ygg,^, to the determinant based 
on the structural completeness. However, This addi¬ 
tional term can be rescaled by a conformal transforma¬ 
tion g^i, —?► (1 + 7)5/ii/ and then can be absorbed into the 
cosmological constant term. In this sense this additional 
term is not expected to affect our results significantly, 
especially at the high energy regime in which the influ¬ 
ence caused by high curvature terms is dominant. In fact, 
one can easily see from the gravitational action that the 
higher order curvatures term will dominate over the ze¬ 
roth order term when curvature gets large. Therefore, 
we will omit this possible additional term in this work. 

Within the Palatini approach we are assuming here, we 
have to vary the action (12.1|) with respect to the metric 
and the connection independently to derive the complete 
field equations. After varying the action with respect to 
we derive the first field equation 

^^q>^''(l + Kl3R)-Xgf^'' 

- ^nPq^Pgo.pg^Pg'^-Rp, = (2.6) 

where q^^y = g^^y + nF^^y and q^’^ is the inverse of q^y. 

stands for the energy momentum tensor. Because 
the matter is assumed to be coupled covariantly to the 
metric g only, the energy momentum tensor is conserved 
like in GR (l^ . 

The second field equation can be obtained by varying 
the action (EH) with respect to the connection 


yf^{aq^'' + Pq^^^ gagg^"") 


= 0 . 


(2.7) 


Note that the covariant derivative Vj, is defined through 
the connection P. 


III. MODIFIED EIBI GRAVITY: A 
PARAMETRIC FRIEDMANN EQUATION 


a-f 4/3 = 1 (2.4) 

to ensure the recovering of Einstein GR at the low cur¬ 
vature limit. Moreover, it can be easily seen that the 


To analyse the behavior of the cosmological solutions 
in the generalized theory defined in Eqs. (Ell and ()2.2|) . 
we follow an approach similar to the one proposed in 
Ref. [13 to our model. More precisely, we will rewrite 
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the field equations in an algebraic form and express the 
quantities of interest using a single variable x, then we 
will represent the behavior of the cosmological solutions 
graphically. We assume q and q~^ denoting and its 
inverse q^’^, respectively. From now on, a hat will denote 
a tensor without making explicit reference to the tensor 
components. We further define (l = g~^q and = 
q~^g to rewrite the first field equation (12.61) as follows 

+ kI 3R) - \{l\^Tr{n-^)Kp{g-^R) - A/ 

= -kT, (3.1) 

where {g~^R) = g^°'Rau{r) and T = T^°‘gav Note that 
I is the identity matrix and TrA denotes the trace of a 
matrix A. 

After taking a trace of both sides of Eq. (13.11) , we get 


Eqs. (ITtI) and (1^ . the energy density as functions of 
LUi and W 2 can be written as: 

4AA^ (^2 — ^i) 

^ q;( 1 + w) — Ar(a;2 — wi)(l — 3i(;) ’ 

where 

q;[q;+ /3(wi + 3a;2)] , „ 

^ —-To- P- 

UJ2 + dwi 


(3.9) 

(3.10) 


Then, we introduce a dimensionless parameter x, 
which satisfies 


Wi = 4 



x 


(3.11) 


|f2|5rr(0-^) = 4A- kT. 


(3.2) 


Therefore, according to Eq. ( 1531 ) . we have 


Moreover, according to the definition of q and Cl and 
Eq. (12.41) we have 


TrCt = 4(1 + Kf]R) + KaR = 4 + kR, 

aK. 


(3.3) 

(3.4) 


Note that equation (13.41) is valid only for a ^ 0. Com¬ 
bining Eqs. dST]), ([32]), dSSl) and (1^ . and after some 
algebra, we obtain 


a{a + iSTrCl) 

Tr{Cl-A 


-(3n+^+ 


l3{a+f3Trrt) 


aX 


7 = 0 , 

(3.5) 


where T = AX/k — T. 

To analyse the solutions within a cosmological scale; 
i.e., we assume the cosmological principle, we first assume 
that the Universe is homogeneous and isotropic at large 
scale and that it is filled with an effective perfect fluid 
with energy density p and pressure p. Then Cl becomes 
a diagonal tensor with 


^0 = Wi, 

n] = uJ2S;. (3.6) 

Therefore, the non-vanishing components of Eq. (13.51) 
read. 


a[a + j3{ijJi -\- 3a;2)]w2 
UJ2 + 3wi 


/3wi 


aX 


/3[a -I- I3{wi + 3a;2)], 


ap 

~T 


(3.7) 


and 


Q;[a-I-/3(a;i-I-3aj2)]wi ap 

-—-/tW2 -I- 

LU2 OCl^i 1 

qA 

= ^ —/3[a-I-/?(wi-I-3a;2)]. (3.8) 


Next, we introduce a constant equation of state w for 
the perfect fluid, i.e. it satisfies p = wp. Combining 



(3.12) 


where Az = {A + KR)/\Cl\i. 

Next, by adding Eq. (13.81) to Eq. (13.71) after being mul¬ 
tiplied by m, the terms involving p are cancelled out, and 
an algebraic equation for x and |f2|3 is obtained: 


Q;(a-I-4/3z|f)|4) . , ImAiI 

1 + 3^4 ^ + -)l^l" 


+I3{a + 4/3z|f2| 4 )(r(; -p 1) = 


aX{w + l)x^ 
(1 -I- 3x^)|fi|3 


(3.13) 


After some rearrangements we derive a quadratic equa¬ 
tion for |0|4 which reads 


Ai|U|5-pi?2|^|3+773 =0, (3.14) 


where 

7?i = Aa/iz(w + x'^) — Piwx^ + —)(1 + 3x^) 

X 

+ Ap‘^z{w + l){l + ^x^), (3.15) 

7?2 = a^{w + x‘^) + a/Aw + 1)(1 + 3x"^), (3.16) 

R^ =-aX{w+ l)x^. (3.17) 


The solution to this quadratic equation, expressed in 
terms of x, reads 




-Ri±\/R2'^-iRiR3 

2Ri 

_A3 
R 2 ’ 


RiA^ 

7?i = 0. 


(3.18) 


In addition, after factoring 7?i(a;) we find that if /3 = 0 
(EiBI limit) or w = 1/3 (radiation domination), 7?i(a;) 
vanishes, which means the second equation in Eq. (13.181) 
is valid. Furthermore, it should be stressed that the ap¬ 
proach mentioned above can not be applied to the case 
in which a = 0 (see Eq. (13.41) 1. Actually, x is fixed to be 
a: = 1 in this case, thus x is no longer a changing variable. 
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However, this fact is not a real problem because one can 
easily derive the cosmological solutions for a pure E? ac¬ 
tion without the need of the approach we are following 

0. 

We have now derived the expression of the energy den¬ 
sity as a function of x in Eq. (13.91) . If we can further 
express the Hubble rate as a function of x, the graphi¬ 
cal relationship between the Hubble rate and the energy 
density can be completed. 

As was already assumed at the beginning of this sec¬ 
tion, we focus on a cosmological symmetry and choose 
a spatially flat Friedmann-Lemaitre-Robertson-Walker 
(FLRW) space-time: 

ds^ = —dt^ + a'^{t){dx^ + dy^ + dz'^), (3.19) 


where t is the cosmic time and a(t) is the scale factor. 
We then define two tensors S and such that 
= V—and E = al -I- l3Tr{Cl~^)Cl. Note 
that h^‘' denotes the inverse of According to the 
field equation ( 123 ), the tensor is the auxiliary met¬ 
ric which is compatible with the physical connection of 

this theory. After some calculations, we obtain |g||E| = 

-' 1/2 

\h\. Inserting it back we have = Sg“^/|E| , and 

The non-vanishing components of 

S and read 


and 


Eg = CTI = a -I- /3(1 -I- 3a;'^), 


E; ^ = 


a + I3{- 


3) 


S" 

3^ 


(3.20) 



hij — fr 1 ^"2CC2 dij d — h‘2SijQj . (3.21) 


Once we have obtained the components of the aux¬ 
iliary metric which is compatible with the physical 
connection of the theory, we can derive the components of 
the connection and the Ricci tensor. After some lengthy 
calculations, we obtain 


Rqo 


^1 p 

h2a^ 


= t( 


3 / /l2 


2 V /i2 


2H] 


3 dho 






(3.22) 


where H = a/a is the Hubble rate and the dot denotes 
the cosmic time derivative. Note that the conservation 
equation p = —3H{p + p) has been applied in the above 
equation. Finally, we arrive at the expression of in 
terms of the variable x: 



a -I- |fi| 4 (4/3z — x^) + 3f^[|H| 4 {x^ — Ajdzx'^) — ax'^] 
3«[2-7|ii/p(1 + ^F 


(3.23) 


In the next section, we will combine Eqs. (13.91) and (13.231) 
together with the definitions (13.6|) . (I3.12|) . (13.201) and 
(13.211) to analyse the past/future asymptotic behavior of 
a FLRW universe in this type of model when filled with 
a perfect fluid with constant w. Eqs. (13.9|) and (13.231) 
correspond to a parametric Friedmann equation being x 
the free parameter. 


IV. THE HUBBLE RATE IN AN EXPANDING 
UNIVERSE 

In the previous section, we derive the expressions of 
p (Eq. (I3.9p l and (Eq. (13.231) 1 as functions of a sin¬ 
gle variable x. Therefore, we can graphically obtain the 
representations of kH^ as a function of np to exhibit the 
behaviors of the cosmological solutions of interest. From 
now on, we will assume a vanishing cosmological constant 
to simplify the analysis, that is, we assume A = 1. 


A. The original EiBI theory: j3 = 0 





FIG. 1. Graphical representation of as a function of Kp 
in the EiBI theory in which /3 = 0. Different curves represent 
different equation of state w. The solid blue, red, green, and 
dashed blue curves show the plots for w = 1/3, 1/5, 0, and 
w = —1/3, respectively. 

As a first glance, we consider the original EiBI theory 
in which /3 = 0 and a = 1. The representations of as 
a function of np are shown in FIG.[TJ One can see that the 
evolution of the energy density terminates at a bounce 
where = 0 and dH/dp ^ 0 at \K,p\ = I for k < 0. 
This bouncing solution is robust against the change of the 
equation of state w. However, if «; > 0 it can be seen that 
the behavior of the Hubble parameter is highly sensitive 
to the choice of w. There are loitering solutions where 
0 and dH/dp — >• 0 at Kp = l/w for w > 0, and 
divergent solutions for w < 0. Furthermore, it can also 
be easily seen that the behaviors of the different curves 
focus around = p/3 when p Ri 0. This property is not 
a surprise because of the prior criteria shown in Eq. (1^ . 
and it can be affirmed in all results shown in the rest of 
this paper; i.e. we recover GR at low energies. Note that 
the results summarized in this subsection are compatible 
with those concluded in the literatures 
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B. Radiation dominated universe 

In this subsection, we analyse if the original loitering 
behaviors and the bouncing solutions within the EiBI 
theory can be altered with the addition of a pure trace 
term to the determinantal Lagrangian, i.e., /3 ^ 0, 
for a radiation dominated universe. The analysis could 
be easily extended to other equation of state but for sim¬ 
plicity we stick to a radiation dominated universe. 


T /3 > 0 





FIG. 2. Graphical representation of kH^ as a function of 
ftp for a radiation dominated universe (w = 1/3). Different 
curves represent different values of p. The solid blue, red, and 
dashed blue curves show the plots for P = 0, 10 “®, and 10 “^, 
respectively. 

We first consider the region in which /3 is slightly larger 
than zero. One should be reminded that in Refs. |37l - K^ 
the authors concluded that the bouncing solutions in the 
EiBI theory for negative ft are robust against the amend¬ 
ment to the EiBI action through an additional f{R) term 
or some functional extents. However, the situations are 
different in our model. One can see from FIG. [5] that 
the bouncing solutions for negative ft are quite sensitive 
to the increase of /3 from zero by even a small amount. 
More precisely, the asymptotic behavior of at large 
p is oc p. This implies the occurrence of a big bang 
singularity in the past. 

On the other hand, for positive ft the loitering effect 
in the EiBI theory becomes a bounce in this model, in 
which oc dp = p— Pmax with pmax being the maximum 
energy density at the bounce. This gives a much regular 
behavior as compared with the asymptotic past behavior 
for ft < 0. 


2. 0 < /3 < 1/4 

As the value of /3 increases and gets closer to /3 ~ 1/4, 
one can see from FIG.[3]that the behaviors of the big bang 
solutions gradually converge to those within the R + R^ 
gravity (/3 = 1/4), that is, = p/3, for negative ft. On 
the other hand, if ft > 0 the loitering behaviors within 


kH^ 



FIG. 3. Graphical representation of kH^ as a function of 
ftp for a radiation dominated universe {w = 1/3). Different 
curves represent different values of p. The solid blue, red, 
dashed blue, dotted blue, and dotted red curves show the 
representations for P = 1/10, 3/25, 7/50, 1/5, and 21/100, 
respectively. The straight line crossing the origin represents 
= p/3, which is the solution within R + gravity 
{P = 1/4) for a radiation dominated universe. Note that 
the vertical lines (from left to right) indicate the maximum 
values of the energy density for P = 1/5 and 21/100, where 
the Hubble rate blows up, respectively. These singular events 
correspond to a big freeze singularity in the past. 

the EiBI theory can be substituted by other interesting 
cosmological solutions. For example, we find that for 
P = 1/10 and P = 3/25, the asymptotic behaviors of H^, 
when p approaches its maximum value Pmax, become 

oc 6p, (4.1) 

where Sp = p — Pmax- Combining it with the conservation 
equation p oc a“^, one can see that this event corresponds 
to a bounce in the past. 

Furthermore, we have also found that the absolute 
value of dH^/dp, which is proportional to H in this model 
as the energy momentum tensor is conserved, is a growing 
function of p. As P approaches P ~ Pi, = 7/50, \dH'^/dp\ 
gets very large at a hnite past cosmic time. Therefore, 
this singular event can be regarded as a quasi-sudden sin¬ 
gularity in the past on the sense that while H is finite, 
H almost blows up in a finite past cosmic time. |f| 

However, if /3 > P^,, i.e. larger than the value cor¬ 
responding to a quasi-sudden singularity, the situation 


We gave in the introduction the definition of the sudden sin¬ 
gularity and the other cosmological singularities related to dark 
energy. In our case this singular event happens in the finite past 
of the Universe. 
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changes drastically. For /? = 1/5 and 21/100, we find 
that the asymptotic behavior of reads 

77^ oc (4.2) 

when p approaches Pmax- Therefore, this event takes 
place at a finite scale factor and a finite cosmic time, 
with both 77 and its cosmic time derivative blowing up. 
These facts highlight the emergence of a finite big freeze 
singularity in the past. 0 

As a summary, we find that the original loitering ef¬ 
fect for positive k can be substituted by a point with a 
minimum scale factor Om, where a bounce (77^ = 0 and 
dH^/dp remains finite), a past quasi-sudden singularity 
(77^ = 0 and dH^/dp nearly diverges) or a past big freeze 
singularity (77^ and dH'^/dp diverge) may emerge. 


3. P <0 





KP 


FIG. 4. Graphical representation of kH^ as a function of 
itp for a radiation dominated universe (w = 1/3). Different 
curves represent different /3. The solid blue, red, and dashed 
blue curves show the representations for /? = —10“^, —1/10, 
and —3/10, respectively. 

In FIG. m we show the representations of kH^ as a 
function oi np for /? < 0. We find that, unlike what 
we concluded previously, the loitering effects (77^ oc Sp^) 
and the bouncing solutions (77^ oc 6p) are robust against 
the decrease of fd below zero. Furthermore, we also find 
that the smaller the value of /3, the smaller the value of 
\kp\ at the loitering event or the bounce. 


I /3>l/4 

On the other hand, in FIG. [S] and FIG. [5] we show the 
representations of k 77^ as a function of k/o for /3 > 1/4 (or 


® Please see Section |I] for the definition of a big freeze singularity 
and the classification of the other cosmological singularities re¬ 
lated to dark energy. In our case this singular event happens at 
a finite past of the Universe. 


a < 0) in a radiation dominated universe. The straight 
line 77^ = p/3, which represents the solution within 7?-|- 
7?^ gravity {jd = 1/4), is also exhibited. One can see that 
the bouncing solutions for negative k are robust against a 
change of fd. However, the loitering solutions for positive 
K in the EiBI theory become big bang solutions where 
77^ oc p for large p. These solutions converge to 77^ = p/3 
when the value of fd approaches f « 1/4. 


k\^- 



FIG. 5. Graphical representation of as a function of tcp 
for a radiation dominated universe (w = 1/3) with f > 1/4. 
The solid, dashed, and dotted curves represent the solutions 
for P = 3/10, 7/20, and 37/100, respectively. Note that the 
solution — p/3, which is the solution within R+R^ gravity 
{P = 1/4), is also shown in this figure. 





FIG. 6. Graphical representation of kH^ as a function of 
up for a radiation dominated universe (w = 1/3) with P > 
1/4 and k being positive. The solid blue, red, and green 
curves represent the solutions for p = 7/10, 4/5, and 9/10, 
respectively. The dashed blue, red and green curves represent 
the solutions for P = 19/20, d = 1 a-nd P = 6/5, respectively. 

Interestingly, we also find from the dashed blue and 
solid green curves in FIG. [5] that for larger values of /3, 
there could be a plateau in the 77^ function for positive 
At, for a radiation dominated universe. This stage may 
correspond to a de Sitter inflationary expansion phase 
after the big bang singularity. This inflationary phase 











is then followed by a classical expansion described well 
in the context of GR. Furthermore, when /3 > 1, the 
solutions with a loitering effect are again recovered (see 
the dashed red and dashed green curves in FIG. (HI)- 
Before concluding, we notice that because this theory 
reduces to GR at the low energy limit, all the radiation 
dominated universe will be asymptotically flat at that 
limit. 


V. CONCLUSIONS 

Since the proposal of the Born-Infeld action for clas¬ 
sical electrodynamics, modified theories of gravity in¬ 
spired on such a proposal and with an elegant determi- 
nantal structure in their actions have been widely in¬ 
vestigated (see Refs. [M3). Despite the large amount 
of works in this subject, the very interesting generaliza¬ 
tion through the addition of a pure trace term into the 
gravitational Lagrangian in the Palatini formalism has 
not been considered before. This modification gives rise 
to the most general action constructed from a rank two 
tensor that contains up to first order terms in the cur¬ 
vature. Such a theory is expected to not only preserve 
the great achievement of GR at low energies, but also to 
generate more drastic deviations from GR than those ac¬ 
complished within the original Born-Infeld-inspired the¬ 
ories at high energies. Modified theories with this term 
have only been investigated in the pure metric formal¬ 
ism [3 and in the teleparallel representation [13. The 
former inevitably suffers from troublesome fourth order 
held equations for the metric or from ghost instabilities 
[3 , which suggests the need of some alternative approach 
to overcome these problems. The latter, which flees from 
the ghosts and results in second order held equations, 
leads in most of the cases to the substitution of the big 
bang by smoother cosmological singularities [13 or a de 
Sitter inhationary stage |15l |. 

Inspired by these motivations, in this paper we gen¬ 
eralize the EiBI theory, which is formulated within the 
Palatini formalism, by adding a pure trace term into the 
determinantal Lagrangian, and analyze the cosmological 
solutions of this theory by assuming a homogeneous and 
isotropic universe for its largest scale. As we expect, the 
early cosmological expansion to be modihed as compared 
with GR or EiBI theory, we assume that the Universe 


is hlled with radiation. Following a similar approach to 
that proposed in Ref. [13 , the behaviors of the cosmolog¬ 
ical solutions are analyzed using a parametric Friedmann 
equation. 

As a summary, we hnd that if k < 0, the big bang 
is substituted by a bounce except for the regions of the 
parameter space 0 < ,5 < 1/4 where the big bang sin¬ 
gularity exists. Note that in Refs. (37l - [39l| the authors 
showed that the bouncing solutions in the EiBI theory 
are robust against the changes of the Lagrangian through 
an additional f(-R) term or some functional extensions. 
On the other hand, if k > 0, we find that the big bang 
singularity can be altered by a loitering effect (/I < 0 or 
/3 > 1), a bounce (0 < ,5 < ,5*), what we named a quasi¬ 
sudden singularity (/3 = /3*), or a big freeze singularity in 
the past < fi < 1/4). However, for 1/4 < /3 < 1, the 
big bang singularity remains. Most interestingly, the big 
bang singularity may be followed by a de Sitter inflation¬ 
ary stage for /3 < 1. This can be verified by the plateau 
in the function as shown in FIG. [i The inflationary 
phase is superseded by a standard cosmological expan¬ 
sion. We summarizes our results in TABLE H] Moreover, 
we should emphasize that the cosmological solutions that 
emerge in this theory are all stemmed from pure geomet¬ 
rical effects. Only a radiation dominated universe is as¬ 
sumed and there is no need of any additional fields or 
exotic matters to drive these cosmological solutions. 
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